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ABSTRACT 



Aims. X-ray observations of gravitationally lensed quasars may allow us to probe the inner structure of the central engine of a quasar. 
Observations of Q2237+0305 (Einstein Cross) in X-rays may be used to constrain the inner structure of the X-ray emitting source. 
Methods. Here we analyze the XMM-Newton observation of the quasar in the gravitational lens system Q2237-I-0305 taken during 
2002. Combined spectra of the four images of the quasar in this system were extracted and modelled with a power-law model. 
Statistical analysis was used to test the variability of the total flux. 

Results. The total X-ray flux from all the images of this quadruple gravitational lens system is 6 • 10"''' erg/(cm--sec) in the range 
0.2-10 keV, showing no significant X-ray spectral variability during almost 42 ks of the observation time. Fitting of the cleaned 
source spectrum yields a photon power-law index of F = 1.82^^55^. The X-ray lightcurves obtained after background subtraction are 
compatible with the hypothesis of a stationary flux from the source. 

Key words, gravitational tensing - X-rays: general - quasars: general 



1. Introduction 

Since its discovery, the gravitational lens system (GLS) 
Q2237+0305 (Einstein Cross; (Huchra et al., 1985) has attracted 
much attention as a unique laboratory to study gravitational lens- 
ing effects. This GLS consists of a quadruply-imaged quasar and 
a lensing galaxy that is the nearest among all known gravitational 
lens systems (the quasar redshift is zq = 1.695, while the lens 
redshift is zc = 0.0395). Due to its very convenient configura- 
tion, it is one of the best investigated GLSs. It has been contin- 
uously monitored by a number of groups from 1992-2005 (Rix 
et al., 1992; Falco et al., 1996; Oestensen et al., 1996; Blanton 
et al., 1998; Bliokh et al., 1999; Nadeau et al., 1999; Wozniak et 
al., 2000; Alcalde et al., 2002; Schmidt et al., 2002). Significant 
microlensing-induced brightness peaks on lightcurves of quasar 
images were detected in this system (see, e.g., Wozniak et al., 
2000; Alcalde et al. 2002; Moreau et al., 2005). The OGLE 
group is continuing observations of the Einstein Cross and the 
database http : / /bulge. princeton.edu/ ~ ogle is being con- 
stantly renewed. 

Observations of the Einstein Cross in X-rays are likely to 
provide important data that can be used to constrain the source 
inner structure. It is hoped that the detection of a variable X-ray 
flux from this system will lead to estimates of the relative time- 
delays At between the images. 

The X-ray emission of Q2237+0305 was first detected dur- 
ing ROSAT observations in 1997 (Wambsganss et al., 1999). 
From the analysis of Wambsganss et al. (1999), a 0.1-2.4 keV 
count rate of 0.006 count/sec was obtained and, assuming a 
r = 1.5 power-law model and Galactic absorption (hydrogen 
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column density Nh = 5.5-10^** cm"^; Dickey & Lockman, 1990), 
a flux of 2.2 ■ 10"'^ erg/(cm^- s) was derived. The 7?(9&4r obser- 
vation of Q2237 -1-0305 did not show any significant flux vari- 
ability, and the spatial resolution of the ROSAT/HRl detector is 
not sufficient to resolve the different quasar images. 

After ROSAT, GLS Q2237+0305 was observed several times 
with the Advanced CCD Imaging spectrometer onboard the 
Chandra X - ray Observatory. Results from a 30 ks and a 10 ks 
observation of Q2237+0305, carried out on September 2000 and 
on December 2001, respectively, were published by Dai et al. 
(2003). For the former of these Chandra observations, the X-ray 
flux was 4.6 ■ lO"'-' erg/(cm^-s) in the energy range 0.4-8.0 keV, 
and the lensed luminosity was 1.0 • 10"*^ erg/s; for the latter ob- 
servation, the flux was 3.7 • 10"'^ erg/(cm^- s) and the lensed 
luminosity 8.3 ■ 10'*^ erg/s (in the same energy range). 

An important point concerns the variability of the X-ray sig- 
nal from the source that carries useful information about the in- 
nermost source structure. In the case of Q2237+0305, the detec- 
tion of X-ray variability is especially important because it may 
lead to a direct estimate of the gravitational time-delays that are 
a significant characteristic of lensed systems. The model pre- 
diction of At between different Einstein Cross images is of the 
order of hours (Schmidt et al. 1998), whereas the optical data 
may provide time-delay values that are accurate to only 1-2 days 
(Vakulik et al., 2006). The only observational estimate of the 
shortest time-delay is obtained by cross-correlating the X-ray 
lightcurves from different images of Q2237+0305 (Dai et al., 
2003) of the Chandra observations: here all four X-ray images 
of this quasar have been resolved and sufficient variability has 
been reported (the difference between total fluxes during the first 
and the second Chandra observations is close to 20% of the av- 
eraged value). The latter observation enabled Dai and collabora- 
tors to determine the time delay of Atba = 2.7;^^" ^ hours between 
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Fig.l. EPIC pn (upper left), MOSl (upper right) and M0S2 (below) 
images of the Q2237+0305. Black circles define areas collecting source 
signal, while white-bordered areas show the regions chosen for esti- 
mates of the background. 



two of the four GLS images. The delays of the other images of 
Einstein Cross are still unknown. 

The XMM-Newton receivers cannot separate different im- 
ages of the Einstein Cross. Nevertheless, one may hope to obtain 
some useful information about the source from the autocorrela- 
tion function of the total X-ray flux from all the images. The 
applicability of this method depends on the presence of suffi- 
cient flux variability and may be complicated by the presence of 
possible microlensing. 

In the present paper we analyze the XMM-Newton observa- 
tion of the GLS Q2237-H0305 taken in 2002. We describe the 
data reduction and present the parameters of the X-ray spec- 
tra obtained by EPIC cameras and RGS spectrometers onboard 
XMM-Newton. Also, we present the results on variability in the 
range 0.2-10 keV. Then we discuss an estimate of the source in- 
homogeneity scale in view of high-magnification microlensing 
events in this GLS. 

2. Data reduction 

GLS Q2237+0305 was observed with XMM-Newton for an ex- 
posure time of 43.9 ks in 28 May 2002 (PI: Michael Watson; 
XMM-Newton pubUc archive). Observation ID, date of observa- 
tion, total source counts and count rates in the extraction regions 
can be found in Table 1 . 

The European Photon Imaging Cameras (EPIC) MOSl, 
MOS2 and pn detectors were operated in full frame mode for the 
whole time of the observation. EPIC data and the data from the 
Reflection Grating Spectrometers (RGS) were processed with 
the XMM-Newton Standard Analysis System (sas, public release 
version 6.5.0).' After filtering of bad-flagged data, the effective 
exposure times were: 39.6 ks for pn-camera, 42.3 ks for MOS- 
cameras and 42.8 ks for RGS-spectrometers. The EPIC maps of 
Q2237+0305 are shown in Fig.l. 

The energy range of EPIC cameras is 0.2-10 keV, and 
the RGS energy range is 0.3-2.8 keV. To obtain the cleaned 
lightcurves, we extracted full count rates from circular areas 
with different radii around the source (see Fig. 1). To estimate 
the backm oLind cfTecls we also used areas of different sizes near 



' SAS software and its description can be found 
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dataset ID 




ni 1 no<^n 1 n t 
Ui lUVDUiUi 


observation date 




Zo.VJ.ZVUZ 






42. J 


exposure tune, (ks) 


pn 


39.6 




I\ ^ TO 


4.9 s 






A no 
U.Uo 






A A/: 


count rate, (count/s) 


pn 






RGSl 


0.11 




RGS 2 


0.10 




MOSl 


3170 




MOS2 


2542 


total source counts 


pn 


11446 




RGSl 


4707 




RGS2 


4279 



Table 1. XMM-Newton observation log for Q2237+0305. Here and 
throughout the paper the name "MOS" is used for the MOS 1 and M0S2 
together. 
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Fig.2. The full band (0.2-10 keV) EPIC pn lightcurves for the source 
(background-subtracted) and the background. The dashed straight Unes 
show the Icr limits of the averaged value of the flux during the observa- 
tion time. The time bin size is 1.8 ks. 
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Fig. 3. The same for EPIC MOSl as in Fig. 2. 
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Q2237+0305 EPIC unfolded spectrum 
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Fig. 4. The same for EPIC MOS2 as in Fig. 2. 

the source on the same CCD plates. Flux variations from dif- 
ferent background areas agree with the corresponding estimated 
Poisson variance. These variations have been taken into ac- 
count in the final error estimates of the resulting net signal. The 
background counts were subtracted (with corresponding factors 
taking into account different background areas) from the total 
counts from the source regions to obtain the cleaned lightcurves. 
The best signal-to-background ratios at low enough levels of 
Poisson variance were obtained for the case of 3-5 arcsec radii 
of the source areas. The background-subtracted lightcurves are 
shown in Figs. 2-4 together with the corresponding background 
lightcurves. 

EPIC spectra were obtained through standard sas proce- 
dures: evselect, arfgen, and rmfgen. For spectral filtering, the 
single high-energy events were extracted in order to identify soft 
proton flares (£ > 10 keV, for pattern zero only; the time-bin 
sizes were 100 s for the pn and 10 s for the MOS). Good time 
intervals were defined using the standard sas procedure tabgti- 
gen (with rate parameter 1.0 cts/s for the pn camera and 0.35 
cts/s for the MOS) and then the event files were corrected. The 
background subtracted EPIC spectra are shown in Fig. 5. 

In the case of the lightcurves, we have not used tabgtigen, 
but we have simply subtracted the background counts from the 
total ones in order to obtain a continuous time series. 

RGS spectra were obtained using standard sas tasks 
(rgsproc and rgsrmfgen). RGS hghtcurves were obtained and 
background- subtracted similarly to the EPIC data (i.e., using an 
evselect procedure); they are characterized by larger errors than 
EPIC lightcurves. No significant absorption or emission lines 
were detected in these spectra. 

3. Spectral fitting 

To fit the Q2237-H0305 spectrum, as a first step we tried a red- 
shifted power-law model (zpowerlw) 

p(E) = KE-^ 

(where F is the photon index), with photoelectric absorption 
(zphabs) by neutral material in the lensing galaxy 




1 2 
channel energy (keV) 



Fig. 5. Unfolded XMM-Newton EPIC (pn and MOS) spectra of 
Q2237+0305 and power-law model. Here the 90% error bars are shown. 
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x'/d.o.f. 


EW, 
keV 


MOSl 


9!-] 




38.8/38 


< 1.6 


MOS2 


8±4 


i.94:;;'|| 


38.6/41 


< 1.1 


pn 


5+^ 


1.85+";[5 


93.0/84 


< 1.3 


MOS 


8±3 


1.92±0.10 


80.1/83 




EPIC 


7±2 


1.88±0.10 


183.5/168 





P(E) 



: Ke 1«G 



Table 2. Absorbed power-law fitting of XMM-Newton Q2237+0305 
spectra - hydrogen column density in the lens galaxy (at a redshift zg = 
0.0395), photon index and I degrees of freedom (d.o.f). Estimated 
EW limits (right-most column, 90% confidence) are calculated for fixed 
values of £ = 5.7 keV and cr = 0.87 keV. The name "MOS" is used 
for the MOSl and MOS2 spectra modelled together; "EPIC" - for both 
MOS and pn spectra modelled together. 



where cr{E) is the photo-electric cross-section and Nh is hydro- 
gen column density (see Table 2), using xspec v. 12.2. 1 . To calcu- 
late lensed luminosities, we used these values for the cosmolog- 
ical parameters: Hubble constant Hq-1Q km/(s- Mpc), Aq = 0.7 
and O-M - 0.3. We compare fluxes and lensed luminosities in dif- 
ferent rest-frame energy ranges obtained from our fit with spec- 
tral fits from Dai et al. (2003) and Wambsganss et al. (1999, see 
Table 3). The row labeled MOS in Table 2 refers to the MOSl 
and MOS2 spectra fitted with the same model simultaneously, 
and the row labelled EPIC represents the MOSl, MOS2 and 
pn spectra fitted with the same model simultaneously. We ob- 
tained the values of Nh for the lensing galaxy taking into ac- 
count the mentioned above frozen value of the Galactic absorp- 
tion (Dickey & Lockman, 1990). The obtained values of photon 
indices and hydrogen column densities are in good agreement 
with the results of Dai et al. (2003), but our power-law indices 
are obviously higher than the power-law index of 1.5, assumed 
by Wambsganss et al. (1999), likely due to the wider spectral 
range. 

In Table 3 we report the XMM-A'ewton fluxes and lensed 
luminosities together with the results from other authors, ob- 
tained from ROSAT and Chandra observations. The fluxes ob- 
tained here are in good agreement with the Chandra ones. 

Commonly considered as evidence of an accretion disk in 
the central engine of quasars or active galactic nuclei (AGN) is 
the detection of a broadened Fe Ka emission line at a rest-frame 
energy of ^ 6.4 keV and a width in the range of 0.1-1.0 keV 
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Satellite 
Instr. 


Energy 
range, keV 


F, 10-'' 
erg/(cm^- s) 


10*5 erg/s 


XMM, MOSl 


0.1-2.4 

0.4-8.0 
0.2-10.0 


2.5+0.3 

5 3+0.4 

-0.3 

6 0+"- ' 


11.8 
13.5 


XMM, MOS2 


0.1-2.4 
0.4-8.0 
0.2-10.0 


26tl 
5.4+0-5 

6 0+" * 


12.9 
14.6 


XMM, pn 


0.1-2.4 
0.4-8.0 
0.2-10.0 


2.4+0.2 

5 l+O 't 

c Q+0.3 

^•*^-n..5 


10.4 
13.3 


ROSAT 


0.1-2.4 


2.2 


4.2 


Chandra 


0.4-8.0 
0.4-8.0 


4.6+"-:^ 


10.0 
8.3 



Table 3. Absorbed power-law fitting of XMM-Newton Q2237+0305 
spectra - fluxes and lensed luminosities. 
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(see Laor, 1991). Such lines were present in the spectra of 42% 
well-exposed objects from a sample of more than 100 AGN ex- 
plored by Guainazzi et al. (2006) and in the spectra of 2/3 of 
the 26 AGN from the sample analyzed by Nandra et al. (2007). 
As shown by Reeves et al. (2004) and Turner et al. (2005), a 
broadened Fe Ka line can also be explained by the presence 
of an ionized absorber with a column density of ^ 10^^ cm"^. 
However, Nandra et al. (2007) show that 1/3 of the cases (in the 
considered sample) cannot be adequately fitted with an ionized 
gas model, and for all of them the line broadening mechanism 
is well explained by reflection from an accretion disk. In our 
case, the apparent lack in Q2237+0305 of prominent absorp- 
tion (mentioned above) suggests that ionized absorber is likely 
not the case for this source. Taking into account the redshift of 
the quasar, the line should be observed between 2 and 3 keV. 
Dai et al. (2003) found the line at 5.7^^ ^^ keV rest- frame energy, 
with 0.87:;:°-^° keV Une width and 1200^^00 gV equivalent width 
(EW), relying on the Chandra (2000-2001) data. 

In the case of the XMM-Newton observation of 
Q2237H-0305, the Fe Ka hue is not evident. We try to 
model this line with the fixed rest-frame energy E and fine 
width a taken from Dai et al. (2003) to obtain an upper limit on 
its EW. The results are shown in Table 2. On the other hand, if 
we let the line energy and the width change over limited energy 
ranges (5-7 keV for the rest- frame energy and 0-1 keV for the 
line width), fitting of the pn spectrum yields E = 6.0;|;j q keV 
rest-frame line energy and a < 0.9 keV for the line width; 
however, this result is not statistically significant (according 
to the test). For MOSl and MOS2 such a treatment does 
not yield a meaningful result at all due to the lower counting 
statistics. 

For completeness, we also fitted the pn spectrum of 
Q2237 -1-0305 with a model often used for the spectral modelling 
of AGN (Soldi et al., 2006, Schurch et al., 2002): a Compton- 
reflected (from neutral material) power-law energy distribution 
with high-energy cut-off Ec (pexrav, Magdziarz & Zdziarski, 
1995): 

p(£) = KE-^e'-k. 

Neutral absorption in the lensing galaxy (zphabs) was also 
included here. The model with reflection provides a slightly bet- 
ter fit than a power-law (x^ = 1.05 vs. 1.18) which is, however, 
not statistically significant. Search for the high-energy cut-ofF 
(Ec) provided no results. The other parameters of this fit are as 



Fig. 6. The distribution of the count differences DY„ = y„+i - Y„, where 
Y„ is the photon counts in every n-th bin. The observed distribution is 
shown with the shaded histogram and the expected distribution arising 
from Poisson noise is presented with the stripped histogram. The error 
bars are at the Icr level. 



follows: photon index F = 2.3+0.3, relative reflection parame- 
ter R = 2.3^ J 2, hydrogen column density in the lensing galaxy 
Nh = 7+^ . 1020 cm-2 and flux of F = 6.0t2l " ^^"^^ erg/(cm2.s) 
in the 0.2-10 keV energy range. 

4. Lightcurves and variability 

The photon counts Y„ in every n-th time bin represent the total 
X-ray flux from all the images of GLS Q2237-I-0305 and have 
been obtained after background subtraction. The time bin size in 
this treatment was varied from T = 0.1 to 2 ks. 

Variability in Y„ (estimated by the ratio of standard devia- 
tion S to the mean Af) is not noticeable in the lightcurves (Figs. 
2, 3, 4). Numerical estimates are presented for 0.5 ks and 1 ks 
bin sizes (Table 4). Note that we do not observe a significant 
correlation (less than 0.2) between pn, MOSl, MOS2 data for 
various bin sizes which indicates that variations of the photon 
counts are simply due to their Poisson dispersion. We examined 
this question in more detail on the basis of the model of station- 
ary Poisson process for the photon count numbers (i.e., with no 
variability). We considered the statistics of the observed count 
differences DY„ - Y„+i - Yn and compared them to the statistics 
of simulated count diflFerences in the case where the differences 
follow Poisson statistics. We find that the distribution DY ob- 
tained from the observational data is very close to our simulated 
DY distribution that assumes a Poisson model for the pn, MOSl 
and MOS2 cameras. The distribution of DY from the combined 
data of all the EPIC cameras is shown in Fig. 6. The;^^^ test yields 
91% probability that the distribution in Fig. 6 is a realization of a 
Poisson process. Also we present in Fig. 7 the structure function 

1 ^ 
1=1 

which agrees with Poisson model. Therefore, the hypothesis of 
constant X-ray flux from the source is compatible with the ex- 
perimental data. The estimates reported above also apply for the 
5-10 keV energy interval, with no substantial changes in vari- 
ability. 
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T, 
ks 


1/ 


y 


pn 

r 


1 


0.13 


0.67 


pn 


0.5 


0.24 


0.89 


MOSl 


1 


0.25 


0.87 


MOSl 


0.5 


0.40 


0.99 


M0S2 


1 


0.33 


1.13 


M0S2 


0.5 


0.40 


0.96 



Table 4. Parameters of the photon counts: T is the time bin size. 
Variability is characterized by V = SIM, S is the standard deviation 
and M is the mean of Y„; the parameter v = S / Vm compares S with 
the calculated Poisson variance. The correlation of the net signal with 
the background in all the cases is less than 0.2. 
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Fig. 7. Structure function (total data from all EPIC cameras). Dashed 
straight Unes show upper and lower Icr-limits estimated for a Poisson 
model. The time bin size is 100 s. 

5. Discussion 

We analysed the XMM-A^ewfo« data (May, 2002) on the to- 
tal X-ray flux from all four images of the GLS Q2237-h0305 
"Einstein Cross". The flux value in the 0.2-10 keV range is 
5.9-10"'^ erg/(cm^- sec), corresponding to a lensed luminosity of 
1 .25 • 10'^'' erg/sec. In order to obtain the true quasar luminosity, 
this value must be corrected, taking into account the lens mag- 
nification. The most recent model estimate of the magnification 
is 16+^ (Schmidt et al., 1998). Fitting the X-ray spectra in the 
0.2-10 keV energy range by an absorbed power-law yields the 
photon index F = 1.88 + 0.10 and hydrogen column density in 
the lensing galaxy Njj - (7±2)-10^'' cm"^. This is in satisfactory 
agreement with the results of Dai et al. (2003). 

The XMM-A'ewton Ughtcurves obtained after the back- 
ground subtraction are compatible with the hypothesis of a con- 
stant flux from the source. An additional argument in favor of no 
variability is the absence of a significant correlation between the 
Ught-curves from different cameras. This does not completely 
rule out source variability: (i) it may be obscured by the Poisson 
variance; (ii) it may be smeared out in the total X-ray flux from 
all the images that have different time delays; (iii) there may be 
periods of different variability amplitudes. These periods may 
correspond to the hfetimes of X-ray emitting blobs orbiting the 
black hole and therefore we might expect variabihty periods of 
the order of the revolution time around the black hole. In any 
case, the net variabihty during the XMM-Newton observation 



cannot be large. As for concern (iii), we note that the ROSAT 
(1997) observations did not reveal significant variability. On the 
other hand, the variability was considerable during the Chandra 
(2000, 2001) observations (Dai et al., 2003). This suggests the 
presence of slow variability changes with characteristic times 
much longer than the inner variability timescale Ty ~ 1 hour 
following from the results of Dai et al. (2003). 

Light-travel arguments imply that the one hour time-scale 
variability corresponds to a size-scale of the X-ray emitting re- 
gion responsible for this variation of Rx ~ TvC/(1 +Zs) = 4 • 10'^ 
cm. An interesting question is whether /?x represents an inhomo- 
geneity in the accretion disk or the entire X-ray emitting region. 
This question might be answered, provided that the black hole 
mass in QSO 2237-H0305 is known. Kochanek (2004) estimated 
this mass, implying a radius of the innermost stable circular orbit 
of 3rg * (4 H- 22.5) • lO" cm (see, e.g., Chandrasekhar, 1983) in 
the case of a non-rotating black hole (where = 2GM/c^ is the 
gravitational radius). The inner radius of a real accretion disk 
must be even larger. However, the mass estimate of Kochanek 
(2004) is model dependent and it is not sufficient to make a 
rigorous judgement concerning a comparison with the order-of- 
magnitude estimate of Rx- 

One may ask how this inhomogeneous structure (if it ex- 
ists) would appear in the presence of microlensing. The theory 
of high-amplification events in the case of a small continuous 
source is well known (see, e.g., Grieger et al., 1988; Shalyapin 
et al., 2002; Popovic et al., 2006 and references in these pa- 
pers). In the case of one emitting spot, the microlensing am- 

— 1/2 

phfication scales as ~ R^^ . Obviously, in the case of a large 
number of small inhomogeneities, the total amplification may 
be smeared out over the source. However, different regions of 
the time-dependent inhomogeneous structure will be differently 
amplified, resulting in an increase of the variability amplitude 
in microlensed image. Note that independent contributions of 
this microlensing effect may reduce the correlation between the 
lightcurves in different images. On the other hand, any promi- 
nent peak in the lightcurve without repetitions after the expected 
time delays can be considered as a probable sign of microlens- 
ing. 

The idea of inhomogeneity of the quasar core in GLS has 
been repeatedly discussed from different viewpoints (see, e.g., 
Wyithe & Turner, 2001; Wyithe & Loeb, 2002; Schechter et al., 
2003). Some evidence in favor of inhomogeneous structure may 
be found in optical observations of different GLSs, where the 
rapid flux variations are observed on the background of slower 
changes (CoUey & Schild, 2003; Paraficz et al., 2006). Note that 
there is an alternative explanation for these variations that invoke 
a planetary-mass object in the lensing galaxy (Schild & VakuUk, 
2003, Colley & Schild, 2003). It might be possible to distinguish 
between these two models by long-term optical and X-ray mon- 
itoring of microlensing events in AGN. 

Observations of the resolved X-ray images with Chandra are 
more informative for the Q2237+0305 variability studies and 
determination of the time delays. Nevertheless, observations of 
the total X-ray flux with XMM-Newton may provide valuable 
complementary information on this issue. In the case of XMM- 
Newton, sufficiently long observations are needed to obtain a 
good result for the autocorrelation function. To avoid uncorre- 
lated contributions from different images to the autocorrelation 
function, the monitoring period must be larger than the differ- 
ential time delays in this system. In the case of Q2237-I-0305, 
the available measurement interval of almost 42 ks is clearly 
not long enough. The model prediction (Schmidt et al., 1998) 
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of maximal delay for Q2237+0305 is 17 hours. Large delays 
(> 35 ks) are also predicted by earlier models of Q2237+0305 
(Rix et al., 1992; Schneider et al., 1988). This is consistent with 
the analysis of the optical lightcurves of Q2237+0305 (Vakulik 
et al., 2006), which yields an upper limit for the delays of about 
~ 2 days. 

While waiting for new observations one may ask whether 
in principle it is possible to determine the variable X-ray signal 
x(t) from the source if the total flux from all the GLS images is 
available. Typically, this problem has an infinite number of so- 
lutions. However, if we know x(t) over a sufficiently large initial 
interval, then x{t) may be recovered uniquely from the total flux 
lightcurve (Zhdanov et al., 2003). This may be of interest if one 
needs to combine results of long GLS monitoring of the total 
flux (cf. XMM-Newton) with shorter observations of separated 
difl'erent X-ray images (cf. Chandra). 
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